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ON : 

We demonstrate analytically that a Bose- Einstein condensate confined in a harmonic trap with 
zero-range attractive interparticle interactions is unstable if there is more than 1 boson. Replacing 
^) , the zero-range interaction by a short-range attractive interaction lifts the instability, and leads to a 

pronounced clustering, by which the particles leak out of the condensate. 
(~^\ \ In a recent preprint jj| , Salasnich claims to confirm analytically the numerical calculations of Edwards and Burnett 
[U and of Dalfovo and Stringari || on the ground state stability of the Bose-Einstein condensate (BEC). Also the 
maximum number of bosons is put forward, for which the BEC with attractive zero-range interaction should be stable. 
However, these conclusions are based on an incomplete analysis of his upper bound K to the ground state energy 
7-H ' (equation (17) in Ref. Q), which reads (in the units h = 1 and m = 1, which are used throughout this paper) 

> : 
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N is the number of bosons, il is the frequency of the confinement potential (assumed isotropic), B is the strength of 
, the zero-range interparticle interaction, and the variational parameter a > is the standard deviation in a Gaussian 
■ trial wave function. 

""^ \ From a mathematical point of view, one can hardly overlook the minimum K — > — oo for a — > in the case B < 0, 
i.e. for an attractive two-particle interaction. The fact that this instability also seems to occur for N = 1 is an artefact 
of an approximation made in the derivation of K. One readily checks that N 2 should be replaced by N (N — 1) . 

The instability with mean distance a — » signifies the collapse of the wave function, and the maximum number 
N max of particles, discussed in Q, should then not be interpreted as the number above which the BEC collapses, but 
as the number of bosons below which a local minimum exists in the Gross-Pitaevskii Q functional. This poses the 
question: how long does it take before the bosons in this metastable state tunnel over the energy maximum into the 
collapsed state, thereby leaking out of the condensate? 

One might object that this singularity in the Gross-Pitaevskii functional in the attractive case is a consequence 
of considering a zero-range interaction, and one might hope that the instability would disappear for a more realistic 
interparticle interaction. However, the Gross-Pitaevskii approach as it stands is hardly applies in this situation. 

Nevertheless, some progress can be made, using the path integral approach which we developed earlier &@|, and 
which allows to determine the thermodynamical quantities, the density and the pair correlation function of a harmonic 
model system. 

By applying the Jensen-Feynman inequality with the trial action resulting from this harmonic model system, one 
finds the following upper bound E v to the ground state energy: 

3 w 2 + il 2 3 (w - il) 2 1 , , f o , , , , 

Ev = jN - { - L + -N(N-1) d 3 vv 2 r )g r 2 
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where w is a variational parameter, v-i (r) is the interparticle interaction potential, and g (r) is the pair correlation 
function of the harmonic model system. The variational parameter w describes the frequency of the 3 (N — 1) internal 
degrees of freedom of the harmonic model system, and takes over the role of the parameter 1 /a 2 in B . 
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In the zero-temperature limit the pair correlation function g (r) as derived in is given by 

,w- (£)"%-*- (3) 

The upper bound E v to the ground state energy can be minimized with respect to w for a rather large class of 
two-particle interactions. In the zero-range case, with v 2 (r, — ry) = — \B\ 6 (r^ — r^), our approach leads to the 
conclusion that E v — > — oo, corresponding to the collapse of the wave function. However, this singularity disappears 
if a finite-range interaction potential is introduced. For instance, a single-step potential of the form 

v 2 (r) = -\V\e(R-r) (4) 

with the range R in the order of a few Bohr radii requires a depth of the order V ~ 10 8 i7 if a scattering length 
a = — 14.5A is assumed Preliminary results reveal that in this case a local minimum is found for w ~ f2 if 
iV < 2875, but that a global minimum for E v is found for w near 10 8 f2, separated from the local minimum by a huge 
energy barrier of order V. 

A fully documented analysis of the (repulsive) Rb case, including anisotropy and temperature dependence will be 
presented in a forthcoming paper ||. For the 7 Li case a detailed study is in preparation. But even without the 
numerical details for the specific alkali metals of interest, the present analysis allows for a few analytical conclusions. 

In the case of a zero-range attractive interaction between the bosons in a harmonic confining potential, the system 
is unstable for N > 2, and the ground state wave function collapses. The minimum in the Gross-Pitaevskii functional, 
found in [jjj, is a local minimum describing the BEC as a metastable state. Our generalization for a finite-range 
interaction, based on the harmonic model system [|5|,[| shows that this instability of the ground state energy is an 
artefact of the 5 function-like interparticle interaction. For more realistic interparticle interactions, this singularity in 
the energy is replaced by a pronounced minimum, corresponding to an extreme narrowing of the wave function, with 
a width of the order of the range of the attractive interaction potential. 
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